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ABSTRACT 

As an alternative approach to the infinite-array description of dense matter in the 
Skyrme model, we report about the properties of a single skyrmion on a compact 
Q \ 3-sphere of finite radius. The density of this matter can be increased by decreasing 

I the hypersphere radius. As in the array calculations one encounters a transition to 

a distinct high density phase characterized by a delocalization in energy and baryon 
charge and by increased symmetries. We will argue that the high density phase has 
to be interpreted as chirally restored one. The arguments are based on the formation 
of complete chiral multiplets and the vanishing of the pionic massless Goldstone 
^ \ modes in the high-density fluctuation spectrum. We show that the restoration of 

I chiral symmetry is common to any chirally invariant extension of the usual Skyrme 

model - whether via higher-order contact terms or via the introduction of stabilizing 
vector mesons which act over a finite range. 



1. Introduction 



^ . In this lecture we will present some studies about the high-density behavior of 

\ baryonic configurations in the Skyrme model II and its variants i' H. These models 

belong to a class of effective models which treat the baryon stabilization (and hence 
the baryon structure) and the interaction between baryons on the same footing. There 
r-| ! is no principal difference between the stabilization and the interaction mechanism, 

the models just act in different topological sectors: in the baryon (=winding) number 
B = 1 sector for the structure physics B and in the B > 1 sector for the interac- 
tion physics i. For B ^ oo one has a new parameter (the density) at ones disposal 
in order to tune the interplay between the stabilization and the interaction mecha- 
nisms. Because of these features the Skyrme model and its variants have the inherent 
possibility to allow for radical changes when baryonic matter is compressed to high 
densities: the baryons which at low densities are well-separated and clearly defined 
objects might completely loose their identity and the baryon matter can become uni- 
form. There are at least two distinct ways for investigating the high-density behavior 
of Skyrme-type models. 

The first approach which was pioneered by Igor Klebanov 11 uses infinite periodic 
arrays of skyrmions with baryon number (winding number) of one per unit cell in real 
space-time. The details of the lattice structure and the periodic boundary conditions 
to which the mesonic fields are subject can be chosen as to avoid nearest neighbor 
frustration at low densities. The aim is to find the classical static field configura- 
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tion which minimizes the energy per cell volume for a fixed baryon charge per array 
cell. The equations of motion (which are non-linear partial differential equations) are 
solved numerically on a lattice grid subject just to two constraints: the choice of the 
crystal and the form of the twisted periodic boundary conditions which guarantee the 
periodic structure and minimize the frustrations in the field gradients (and thus the 
energy) between neighboring cells in the asymptotic low density region. The bound- 
ary conditions are then extrapolated without alteration to high densities as well. At 
low densities the skyrmions in the periodic arrays are well-separated, ensembled in a 
phase of weakly interacting baryons. With increased density, however, they grow in 
size as measured by their r.m.s. radius until at a critical density - as first observed 
by Wiist, Brown and Jackson - they "melt" to a distinct high-density phase where 
the skyrmions completely lose their identity. The new phase is characterized by the 
approximate uniformity in the baryon as well as the energy density, by the fact that 
the averages of the a- and vr-fields vanish over the cell volume and that there appears 
an additional symmetry at the critical density: the "half-skyrmion symmetry" iof 
Goldhaber and Manton.ll 

The short-comings of these calculations are: they are by fiat of numerical nature, 
the rotational and translational symmetries are explicitly broken down to discrete ones 
by the artificial crystalline grid structure. Finally the calculations depend on the form 
of the crystal and the twisted periodic boundary conditions which are extracted in 
the asymptotic low density regime, but nevertheless applied at high densities. For 
instance, the order of the phase transition can depend on these choices: For the 
same simple cubic lattice it is first order for Klebanov's choice of boundary terms, 
while it is second order for the modified rectangular boundary terms of Jackson and 
Verbaarschot H where a preferred direction in space is singled out. But the fact that 
there are phase transitions are common for all these calculations - and even the critical 
densities are approximately the same. Klebanov's crystal is not of minimal energy, a 
face-centered cubic arrangement of skyrmions (at low densities) is preferable. Again, 
as the density increases there is a second-order phase transition and the half-skyrmion 
symmetry emerges.liiJ The crystal of minimal energy is in this higher density phase 
and has energy per baryon only of 3.8% above the topological lower bound. Further 
generalizations are the inclusion of temperature and the derivation of the equation 
of state for the Skvrme matter in the case of the simple cubic lattice 1111 as well as 
the fee structure by T.S. Walhout. He has also studied periodic arrays using an 
a;-stabilized variant of the Skyrme model.lll 

Nevertheless, the question might arise whether these phase transitions are lat- 
tice artifacts or not. Fortunately, there is an alternative approach to study the 
same phenomena by replacing the periodic arrays in fiat space by few-skyrmion- 
systems (or even a single skyrmion) on the compact manifold S^(L) pioneered by Nick 
Manton. The finite baryon number on S^(L) corresponds to a finite baryon den- 
sity in R^(with infinite baryon number), so that one obtains a model for skyrmion 



"This symmetry is characterized by famihes of planes on which the scalar field, a is zero. The pionic 
fields have a reflection symmetry about these planes while the scalar field is reflection antisymmetric. 
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matter which is appreciably simpler to study than any lattice model. The density of 
this matter can be increased by decreasing the hypersphere radius L, the radius of 
the 3-sphere in four dimensions. In the limit L — oo a skyrmion localized on the 
hypersphere has basically the same properties as an isolated skyrmion in flat space 
or in a periodic array for infinite separation. For a large hypersphere it is still well 
localized and dominated by its stabilization mechanism, but because of the curvature 
and finite size effects of the hypersphere it acts as if the tail of another skyrmion (here 
of course nothing else than itself seen via the opposite pole of the 3-sphere) is present. 
By shrinking the hypersphere radius this interaction mechanism will become stronger 
and stronger. Finally it can become so strong that the stabilizing mechanism cannot 
localize the skyrmion any longer: the skyrmion will be smeared (=delocalized) over 
the whole 3-sphere. The connection between both approaches can be made by iden- 
tifying the averaged baryon densities: in the periodic array calculation the baryon 
density is averaged over the cell volume, on the hypersphere it is given by the ratio 
5/(27r^Ljjij^) of the baryon charge B and the (surface-) volume of the 3-sphere in 
terms of the dimensioned radius Ldim- The S^(L) approach has the advantage that it 
is technically far simpler than the period array one Band that it allows for mathemat- 
ically rigorous results for the properties of the "high- density" (or "small L") phase 
of a single skyrmion on S^(L). Furthermore neither the (continuous) rotational nor 
the translational symmetry are broken. The disadvantage is obvious: The physical 
space is not a hypersphere - at least not a small one. Both approaches supplement 
each other in a complementary way: the first is realized in normal space-time and for 
infinite systems (a precondition for a phase transition), but is technically complicated 
and plagued by its crystalline nature, by the ambiguities in choosing the array itself 
and the boundary terms; the other is simple and exact, but is set up in an unphysical 
world and limited to finite systems which therefore should be not identified with an 
infinite periodic array, but just with one array cell. 

The talk is organized as follows. In sect. 2 we will give Manton's general description 
of the Skyrme model on the 3-sphere and review some of the B = 1 properties. In 
sect. 3 we present arguments why the phase transition has to be interpreted as chiral 
restoration. In sect. 4 we report on the generalization to variants of the Skyrme- model 
which involve higher-order contact terms or which are stabilized by vector mesons. 
Sect. 5 contains a discussion about the order of the phase transition. We end the talk 
with a short discussion section. 

2. The Skyrme Model on the Hypersphere 

As mentioned there is considerable interest in studying the Skyrme model 



on a 3-dimensional hypersphere of radius L, S^(L).t4'0JI We will report here espe- 

^In fact most of the calculations can be done analytically or involve at most the numerical task of 
solving ordinary non-linear differential equations. 
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cially on the findings of refs. In order to get parameterization-independent 

results it is useful to relate the quaternion representation U of the Skyrme model 
to a cartesian representation = (<l'°, "l*^, with the imposed constraint 

= 1 as f/ = $° + ir ■ $ and to introduce the strain tensor HI 

Kij = di^'^dj^'' (with = ((T,7r^7r^,7r2^)) 

= -\TT{Wd,U,Wd,U}+, 

where i and j label the space coordinates. This is a symmetric 3x3 matrix with 
three positive eigenvalues which we will denote as A^, Af, A^. Manton has shown 
that the Aj have a simple geometrical interpretation. They correspond to the length 
changes of the images of any orthonormal system in a given space manifold (here 
or S^(L)) under the conformal map, U, onto the group manifold, here SU{2) = 
S3. Therefore the name strain tensor which refers to such a general "rubber-sheet" 
geometry. Among the invariants of Kij, three are fundamental Hand have a simple 
geometric meaning: 

Tr{K) = XI + XI + XI = E length' 

|{(Tr(ir))'-Tr(ir')} = A^Ag + A^A^ + A^A^ = Earea' (3) 
det{K) = XlXlXl = volume' . 

The first one measures the sum of the squared length changes of the mapped or- 
thonormal frame, the second one the sum of the squared area changes and the 
third the squared volume change. With the help of Eq.(0) it is easy to see that 
these invariants are (modulo normalization factors) the static energy densities of the 
second-order non-linear a model, (/'/2)Tr(9jf/(9jf/'^), the fourth-order Skyrme term, 
— {el/4)TT[WdiU, WdjU], and the sixth-order term proportional to the square of the 
baryon density (see Eq. (p4D in section 4). The static energy of a general skyrmion 
configuration on the hypersphere has in this language the form i 



E = f,^^dV{Xl + Xl + XI + xlxl + xlxl + xlxl) 

= [ dV (Xa - XbXc? + (cycl. perm.'s) + 6 / dV XaXbXc- (4) 

Js'^iL) JS3{L) 

The last term is just 127T^B with the baryon (winding) number B, since AaAbAc/(127r') 
is the Jacobian of the map S^(L)— s>S^(l)= SU{2) - in other words the baryon (wind- 
ing) number density. Because of the positive definiteness of the terms in the last 
expression of (^ it is obvious that any skyrmion configuration has to respect the 
topological bound 

E > 127r'5. (5) 



'^All other invariants can be constructed from them. 

''We have adopted a dimensionless form of the Skyrme lagrangian. To obtain the dimensioned 
quantities, one divides the A^'s and multiplies L by 2^/264/ f.^ and multiplies E by ^f^t^j.^. 
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Furthermore, one can immediately see that there is exactly one possibility to satisfy 
the topological bound: The insertion of the identity map, = A;, = Ac, with Aj = 1 
which corresponds to an isometric mapping from the spacial manifold into the target 
manifold and which has baryon (winding) number B = 1. Specifying to the case at 
hand where the target manifold is SU{2) = S^(l), the unit 3-sphere, we see that 
the only way for saturating the topological bound (the absolute minimum of any 
skyrmion configuration with a non-zero winding number) is the isometric mapping of 
the spatial S^(L) with L = 1 onto the target iso-spin sphere S^(l). We can therefore 
conclude that the topological bound can nether be saturated by any B > 1 system or 
by any skyrmion configuration in fiat space, R^.iAfter minimizing E[Xa, Xt, Ac] under 
the constraint 

dV XaXbXc = 27i'^B (6) 

53(L) 



Manton li^ and Loss Il3 found that for L < 1 there is even a stronger bound on any 



B = 1 skyrmion configuration, 

^identity = OvT^ + ^) , (7) 

which is still satisfied by the identity map Xa = Xb = Xc = I / L independently of the 
parameterization of the B = 1 configuration. In other words for L < 1 the identity 
map is the absolute minimum of any B = 1 configuration. 

Finally, it is now simple to show in full generality that for L < ^2 the 

identity map is stable for all allowed {6B = 0) small-amplitude perturbations Aj —>■ 
(1/L) +6i{x), whereas it becomes a saddle for L > V2.0[Jnder the small perturbation 
given above the energy is given to second-order as 

E[Xa, Xb, Ac] = ^identity + + Js) j + + J2) j + J2 j + ^i^') 

with Ii = 6a + Sb + 6c, I2 = 6l + 6"^ + 6l and J3 = 6a6b + + 45a- Using the 
5-constraint (^ we have up to third order corrections the relation 

^/ dVh = -\( dVh + 0{6^). 

After inserting this relation in Eq.(|^) and using the fact that the 5j's are local vari- 
ations, one can finally find that the identity map is stable against small-amplitude 
fiuctuations, i.e. E[Xa, Xb, Xc] > -^identity) as long as the following local inequality 
holds 

1 + J2) + + ^l) - 2('^a5fe + 6b6c + 6c6a) > 0, (9) 



"^Of course, we cannot exclude the possibility that such a configuration might come infinitesimally 
close to the topological bound. 

■'^That means that the identity map is at least a local minimum for 1 < L < V2. 
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in other words as long as L < This concludes Manton's general proof - indepen- 
dent of any parameterization - that the identity map is at least a local minimum up 
to Lc = y/2 and becomes a saddle for L > a/2. 

So far we have not specified any parameterization. In the following we will spe- 
cialize to the familiar hedgehog parameterization. It has the following properties: 
(a) according to a general theorem by Palais llj about reduced variational equations 
for symmetric fields a variational solution of hedgehog form is a solution of the full 
set of Euler-Lagrange equations of the Skyrme model, (b) the hedgehog solution is 
the energetically lowest solution known so far in the B = 1 sector in and (c) 
it is numerically proven that the B = 1 hedgehog solution is at least a local min- 
imum. Hedgehog configurations on the hypersphere are conveniently described in 
terms of the conventional 'polar' coordinates fi, 9, (p on S^(L) with < /i, 6* < vr and 
< < 27r, so that a typical point on S^(L) has Cartesian coordinates 

(Lcos/x, LsiuficosB, L sin fi sin 6 cos (p, L sin fi sin 6 sin (p). (10) 

In these coordinates the metric is ds"^ = L'^i^djjp' + sin? jjidO^ + sin^/i sin? dcpP') and 
the volume element is 

dV = sin^ jjidn sin dO d(j) . (11) 

The field components of a hedgehog of baryon number B on S^(L) have then the 
form 

$0 = cos/(/i) 

= sin /(/i) cos ^ 

= sin f{fi) sin 9 cos (p 
<l>^ = sin /(/i) sin 61 sin , (12) 

where the "radial" profile function /(/i) is subject to the boundary condition /(O) = 
and /(vr) = Svr, B integer. The topological winding- number B is of course the baryon 
charge of the corresponding field configuration. The field components, /tt^" (/tt is 
the pion-decay constant) should be identified with the familiar cr, vr^, vr^, vr^ fields. 
In the normal quaternion representation the hedgehog ansatz reads 

U = ^^ + iT-$ = exp(zf • f{e, (p) fifi)) (13) 

where f(6', (p) is the usual radial unit vector. Note that for the hedgehog con- 
figuration the eigenvalues of the strain tensor are simply = /'^/L^, = = 
sin^ //(L^ sin^ /i) where /' stands for 

The energy of the hedgehog field configuration ([13D on S^(L) is li 
E = AttL r sin^ /irf/i { f + 2-''''' 

JO \ J 



sm^ /i 



1 /""' 

-|- A-K— I siv? jjidjji 
L Jo 



sin^ fi 



2r- + ^ 

sm /i 



(14) 



^Note in order to obtain the dimensioned quantities, one multiplies L in ( |T^ ) by 2\f2e^l j.^ and E 
by V2e4U- 
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Variation of (JT^) with respect to the radial profile function /(//) leads us to the Euler 
equation 

d^/ mi2ndf sin 2/ 
d/x^ sin^ yU sin^ /i 



sin /i rf/i^ sin /i 1 \ ci/i/ sin /i 



Again because of Palais' general theorem solutions of ( ]T5| ) yield solutions of the 
full set of Euler-Lagrange equations of the Skyrme model on S^(L) when substituted 
into (0) or (|r^). It is very simple to check that one solution of Eq.(]TB|) (in the 
case 5 = 1) is the map /(/u) = ^ which corresponds to the uniform mapping of 
S^(L) onto the isospin manifold S^(l) = SU{2). One can easily see that for the 
uniform hedgehog map, /(/i) = /i, the eigenvalues of the strain tensor (0) are simply 
= = = In other words it is a special parameterization of the identity 

map. The uniform map, /(/u) = /i, is a solution for all L, but has no finite L ^ oo 
limit. The energy associated with the uniform map is of course given by (|^). As 
mentioned for 1 < L < -\/2, the uniform map is still at least a local minimum. For 
L > \/2 there exist a lower energy configuration which represents a skyrmion localized 
about one point, i.e. f{fi) 7^ in the hedgehog parameterization In fact 

in this parameterization the identity map bifurcates into two solutions of the same 
energy concentrated around the north or south pole of the hypersphere which are 
related by 
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/jv(/i) = 7r-/5(7r-/i). (16) 

In the limit L 00, one recovers the usual fiat space skyrmion.0 The identity 
map has 0(4) symmetry {see Eq.(|T2p) and is therefore completely uniform in energy 
and baryon density. For L > \pl the symmetry of the energetically most favorable 
solution is broken to 0(3) indicating that the solution is localized in energy and 
baryon distribution. The bifurcation at L = a/2 is of the standard "pitchfork" type 
corresponding to a second order phase transition. In Landau's scheme the bifurcation 
is characterized by a symmetric quartic polynomial whose quadratic term changes 
sign. 

In order to discriminate the delocalized (high-density) phase from the localized 
(low-density) phase in a way which can be generalized to few-skyrmion systems on 
the hypersphere or to flat space array-calculations, several candidates for an order 
parameter were studied in ref. 0: One candidate is the integrated squared deviation of 
the local energy/baryon density from the averaged density normalized by the squared 
averaged density times the hypersphere volume. It converges to the value 1 for the 
localized solution in the limit L — » 00, for smaller L it decreases and shows a square- 
root bifurcation at L = -\/2 to the value zero belonging to the identity map. The 
interpretation would be that this order parameter signalled deconflnement in the 
energy and baryon distribution. This is, however, a special property of the B = 1 
system on S^(L). For few-skyrmion systems 011 and flat space arrays this 
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parameter never becomes zero in the high-density phase, although there is still a 
sizable decrease in its value from the low-density to the high-density phase. The 



second suggestion of ref. 113 on the other hand works also in these cases. It is the 
squared chiral expectation value 

{<y/u? + (v?/^)^ (17) 

where the a and pion fields are averaged over the hypersphere volume (or the cell 
volume for periodic arrays). Because of the residual 0(3) symmetry for the hedgehog 
(or discrete symmetries in the array calculations) which cause (tt) to be zero it has 
the same content as the parameter (a) alone. For the identity map it is obvious 
(see Eg . ([I^) ) that the parameter ([T7|) vanishes. In the case of the localized solution 
it approaches the value 1 with increasing L, since in most of the space the U field 
of the localized skyrmion is approximately close to the vacuum value Uq = 1 and 
only markedly deviates from this in the small region where the skyrmion is located. 
The interpretation is that a zero in this order parameter signals chiral symmetry 
restoration or at least "chiral democracy" , since there is no escape from the "chiral- 
circle" constraint cr^J- vf^ = But nevertheless as shown for few-skyrmion systems 
on the hypersphere 0' and for fiat space arrays S' the a and pion profiles of 
the high-density phase are so equally distributed about the chiral circle that their 
averages vanish. 

Let me summarize the B = 1 results on S^(L): There is a second-order transition 
at a baryon density l/(27r^Lc^) for Lf. = V2 from a localized and chirally broken low- 
density solution to a "chirally restored" high-density solution. Restoring dimensionful 
parameters into the Skyrme lagrangian (by using the empirical value of the pion decay 
constant = 93 MeV and an = 0.0743 guaranteeing a reasonable value for qa) 
one finds 



^ ^ 0.20fm" 



\2V2eJ v^'27r2 

which is far too low .0( remember that the nuclear matter density is /Onm = 0.16fm~^). 
The important result, however, is that the S^(L) value for p^. is more or less the 
same than the one found in the periodic array calculations under the same 
input parameters, approximately pc = 0.17fm~^. Furthermore, the corresponding 
transition densities pc in the few-skyrmion calculation 0' HI interpolate between these 
two numbers. Thus the identification of the hypersphere formalism with the periodic 
array calculation seems to be justified even quantitatively. 



'This might improve if the stabilization term is replaced by more realistic ones. Furthermore all 
contributions from kinetic terms are neglected. Finally, there is the problem of the missing central 
attraction ^ which is probably linked to loop corrections ^ or missing 1 /Nc corrections ^ . Any 
additional attractive term tends to localize the skyrmion and will therefore move the transition point 
up to higher densities. 
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3. Indications for a Chiral Symmetry Restoration 

In this section I will present some arguments (see refs. EH El for more details) 
why the transition reviewed in the preceeding section should be interpreted as chiral 
symmetry restoration. As mentioned the complete uniformity of the high-density 
B = 1 solution on S^(L) is linked to the strong condition /(/x) = and does not 
generalize to few-skyrmion systems 0' HI on S^(L) or periodic arrays ME in fiat 
space. However, the "chiral democracy" as discussed at the end of the preceeding 
section can even be achieved via a weaker condition 

/(/i) = 7r-/(7r-/i) (19) 

which indicates that in the high-density phase there is a symmetry about the hyper- 
sphere equator between the "northern" and the "southern" hemisphere. In fact this 
is just the half-skyrmion symmetry H which is a common signal of the high-density 



phase for few-skyrmion systems on the hypersphere lllr 11^ as well as for periodic ar- 
rays in fiat space i' H. So the chiral symmetry restoration seems to be prior to 
the delocalization. In fact this can be tested by adding a (pion mass) term which 
explicitly breaks chiral symmetry:!! 

/:sB = -^Tt{U + W- 2). (20) 

When this term is added, the phase transition is not a sharp one any longer, but 
it is smeared out. The chiral order parameter (|T^ still approaches zero for higher 
and higher densities, but never actually becomes exactly zero. Of course, this is just 
a consequence of the explicit breaking of the chiral symmetry. A second effect of 
adding the term (|20|) is a small shift of the now approximate phase transition to 
higher densities, since the term induces an attractive force. 

In ref. H it is furthermore argued that the vanishing of the order parameter (a) 
(for the Skyrme model without explicit symmetry breaking of course) implies the 
vanishing of the matrix element for pion decay {0\A1^{x)\7i^) in case the axial current 
is calculated in the mean-field approximation in the frame-work of a cx-model, 
and therefore the vanishing of the "effective" pion decay constant and the quark 
condensate. 

The most convincing arguments why the phase transitions in the Skyrme model 
ought to be identified with a chiral symmetry restoration are the following two: In 
the B = 1 case on S^(L) it can be shown that the transition from the low-density 
phase to the high-density phase is accompanied 

(i) by the formation of complete chiral multiplets in the fluctuation spectrum 

(ii) and by the vanishing of the three pionic Goldstone modes from the spectrum. 

These two points are treated at length in ref. H where the local stability of the 
uniform anc? the localized skyrmion solution on S^(L) and thus their small-amplitude 
normal modes are investigated. If all normal modes have positive energy, the solution 
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is locally stable. If some normal modes have negative energy, the solution is a saddle 
point. The study of the small-amplitude fluctuations about the B = 1 solutions on 
the hypersphere revealed the following scenario: 

(i) In general, the modes about a hedgehog solution have just the usual 0(3) de- 
generacies, i.e. a degeneracy 2N + 1 in terms of the principal quantum number 

> 1 of the modes. 

(ii) The modes about the identity map on the other hand can be classified according 
to the following representations of the group 0(4): the symmetric tensor repre- 
sentation (A^, 0) which has (A^ + 1)^ degenerate states and as static eigenvalues f 

A^,o = L {N{N + 2) - 4} + 2 {N{N + 2) - 2} /L, (21) 

and the (A^, 1) representation with a 2((A^ + 1)^ — 1) degeneracy and the static 
eigenvalues 

Xn,i = {N{N + 2)-3}{L + 1/L). (22) 

Thus their degeneracy exceeds by far the one of the modes about locahzed 
hedgehog solutions. Furthermore for each mode about the identity map with a 
given parity, there exists at least one degenerate mode of opposite parity. Using 
the covering group SU{2)l x SU{2)ji of 5*0(4) it can be shown (see ref. c3) 
that the (A^, 0) modes belong to a trajectory based on the four-fold degenerate 
(1/2, 1/2) representations of SU(2)l x SU{2)r, while the (A^, 1) modes belong 
to a trajectory based on the six-fold degenerate (1, 0) -|- (0, 1) representation. 
Thus the modes form complete multiplets of chiral symmetry either by the 
additional (cr, vf) degeneracy of the (A^, 0) modes (|2T| ) or by the parity doubling 
of the (A^, 1) modes (P^. Note that in ref. El the (A^, 0) modes were identified as 
purely "electric" grand-spin modes (but with a degeneracy between e.g. grand 
spin K = O"*" and K = 1~ modes), while the (A^, 1) modes were identified as 
degenerate "magnetic" and "electric" modes of the same grand spin K, but 
opposite parity. The (A^ = 1, 1) modes are the sixth zero modes about the 
identity map. 

(iii) At high densities (L < \/2) all normal modes about the identity map have 
positive energies, the identity map is stable. There are four degenerate low- 
lying modes in the spectrum, the four (A^ = 1, 0) modes, one with a quantum 
numbers and three with pion ones indicating that the a and all three vr-fields 
are degenerate and treated on the same footing. They have the static eigenvalue 

A;v=i,o = -L + 2/L. (23) 

(iv) The four (A^ = 1, 0) modes cross zero at L = a/2 and become negative for 
L > \f2 signalling there the instability of the identity map. 

*We use the phrase "static eigenvahie" as a short-hand notation for the change of the static energy 
under the smaU perturbation of the fluctuation-niode. 
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(v) At L = \pl there is the bifurcation into the normal locahzed hedgehog solution 
with the reduced 0(3) symmetry. All modes about the localized solution are 
stable for L > a/2. The modes with N > 1 have (with the exception of some 
accidental degeneracies for special values of L) only the usual 0(3) degeneracy 
2A^ + 1. However, there exist 9 zero modes.l3The three additional zero modes 
(compared with the number of zero modes for the identity map) correspond to 
the three unstable pionic {N = 1, 0) modes of the identity map. Their fourth 
partner, the cx-like mode, corresponds to a positive stable excitation about the 
localized hedgehog (to an infinitesimally small conformal variation in the profile 
function f{fi) /(/i) + The breakdown of the continuous symmetry 
0(4) of the identity map at L = to the 0(3) symmetry of the localized for 
L > \/2 energetically favorable hedgehog solution is therefore closely linked to 
the appearance of three new zero modes. 

(vi) These three zero modes should be in fact interpreted as the three pionic Gold- 
stone modes linked to the spontaneous breakdown of chiral symmetry. The 
arguments for this identification are the following: Whereas six of the nine zero 
modes stay normalizable even in the limit L — > oo, three have no finite L —>■ oo 
limit. This is consistent with the fact that the fiat space hedgehog has only 
six zero modes. Furthermore when the chiral symmetry breaking term ([2^) is 
added to the Skyrme lagrangian, the above mentioned normalizable zero modes 
stay zero modes, whereas the three additional zero modes are shifted in energy 
by the pion mass. These are the three pionic Goldstone modes which become 
massive when an explicitly chiral symmetry breaking term is added and which 
show up as non-normalizable plain wave excitations of the vacuum in the flat 
space limit. 

Since the three Goldstone modes disappear for the high-density delocalized phase 
and the high-density modes form complete chiral multiplets, the chiral symmetry 
restoration is established. 

4. Generalizations 



4.1. Contact Terms 

In the following we will show that the form of the Skyrme model (|l|) and especially 
the stabilization by the fourth-order Skyrme term is not a necessary precondition 
for the very existence of the chiral phase transition. Let us for instance replace 
the fourth-order stabilizing term in (|lD by a stabilizing term of sixth-order {see e.g. 
ref. 1^3) 

£2,6 = ^Tr [d.Wd'^u) - c,B^B^ (24) 



-'Note that a hedgehog in flat space has only 6 zero modes: 3 rotational (which cannot be separated 
from iso-rotational ones because of the hedgehog symmetry) and 3 translation ones. 



11 



(note Cq > 0) where is the topological baryon current 

= ^Tr {{WdM) (U^d^U) {WdpU)) . (25) 

As reported in section 2 the fourth-order Skyrme term as well as the sixth-order term 
and the second-order non-linear sigma model term have a geometrical meaning when 
expressed in the strain tensor language. One can apply therefore Manton's general 
machinery to show that the sixth-order term allows for the identity map as a 
solution and to construct the critical value where the identity map of the model 
( |2^ ) becomes unstable and bifurcates into a localized solution. In fact all qualitative 
phenomena show up as before. Naturally the value of Lmin for the minimum of the 
static energy is in general not the same for both models and more importantly the 
critical hypersphere radius where the chiral restoration occurs is 

Le = (3)3L^i, (26) 

for the model ( p4D instead of = \/2Lyam as it was the case for the usual Skyrme 
model (|1|).0 



As shown in ref. c3 any lagrangian which (a) is a polynomial of the fundamental 
invariants (^ of the strain tensor, (b) which has the usual vacuum properties and 
(c) which allows for (locally) stable B = 1 Skyrmion solutions, leads to a chirally 
restored phase at high densities. In addition the position of the phase transition can 
be uniquely derived just from the lagrangian without solving any equation of motion. 
Let us discuss for this purpose a lagrangian of the analytic form 

oo 

Cg ^ -g{-C2,C,,C,) = - J2 ^nnh {- {- {- C.Y' (27) 

ii, i2, j3=o 

where ji,j2,j3 are integer indices and Aj-^j^j^ Taylor coefficients. The arguments of 
the analytic function Q and the signs are chosen in such a way, that the corresponding 
static energy density has the simple form 

oo 

£ = Qi£2, £a, £6) = E "^nhjs £2''£^£^- (28) 

The function Q should be of course positive definite to ensure a positive energy density 
and ^(0, 0, 0) = {e.g., Aqqq = 0) to ensure triviality in the B = sector. By inserting 
/ = or /' = one furthermore learns that Q{x, 0, 0) > and Q{x, 0) > for 
any real x > 0. The first inequality signals that a possible negative term expressed 
solely by the second-order term can nether be compensated by adding fourth- or 
sixth-order terms in any combination to the static energy density. Therefore Skyrme- 
type models which at fourth-order are not positive definite cannot be saved by the 
addition of a sixth-order term proportional to B^B'^ or any power of this term. 

Note that the energy density ( pH] ) is the most general symmetric function of the 
eigenvalues, A^, of the strain tensor Kij and includes all forms which can come 
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form lagrangians which involve arbitrary combinations of powers of first derivatives 
of the fields to even order. 

In the chirally restored (uniform) regime on the hypersphere the energy densities 
of the second-, fourth- and sixth-order term, £2, ^4 and Sq, have the simple form 
(modulo a prefactor which can be incorporated into the Taylor coefficients ^(Jij2j3)) 

£2^ = ^, i = l,2,3. (29) 



In ref. the eigenvalues of all possible small amplitude (static) perturbations around 
the identity map were obtained as the following analytical result: 



- / V A ^ +A 



with 



KN,a = I [n{N + 2) + - 4 - ^ ({a(iV + 1)}^ - 4) } L'-'' 




d 



882 L^dSi L^dS, 



i - 1 



where the indices > 1 and a = 0, 1 characterize the allowed 0(4) representationsiH 
{N,a = 0), the symmetric tensor representations which have (A^ + 1)^ degenerate 
states and the {N,a = 1) representations which have a 2{(A^ + 1)^ — 1} degeneracy 
(see refs. HI H for further details). Note that in the case a = 1 the expression for 
X%^a=i simplifies to 

A^,,=i = L{N{N + 2) - 3} J- + -1 J-^ g{S2,S„S,). (32) 

A negative value for any eigenvalue Aat^i in (P^) would mean that infinitely many of 



the (iV, 1) modes would pass through zero energy simultaneously. This is probably 
indicative of a non-perturbative configuration of lower dimension, e.g., a string-like 
configuration. Fortunately, one can exclude this pathological instability, since it vio- 
lates the constraint that the model (|27|) should have the same vacuum properties as 
the nonlinear a model (see ref.ll^ for a discussion about this point). So the (A^, 1) fiuc- 
tuations cannot lead to an instability. In case the identity map is (locally) stable in 
the high-density region {i.e. all A^^ > 0), the monotonic increase of the eigenvalues 
of the small-amplitude normal modes with increasing is guaranteed. Then, 
the critical Lc where the identity map becomes unstable against small- amplitude per- 
turbations is given by the value of L where the first A^q in Eq.(pOD becomes negative. 
Finally, let 

g{S2,S^,£e) for^2i^0 (33) 

and 

g{S2, ^4, Se) ^ CiS2^ + C2^4 for S2^ ^ oo (34) 



See section 3 for the corresponding results for the usual Skyrme model ffl) 
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where i runs from 1 to 3 and Ci > and C2 > are fixed positive coefficients. The 
condition (^) enforces that asymptotically for low densities (large L) the lagrangian 
( P7| ) is dominated by the non-linear sigma model lagrangian, whereas for high densities 
(small L) it scales as a free Fermi gas {i.e. Sg oc 1/L^). The low density behavior 
is of course indisputable, the high density constraint on the other hand requires the 
additional input that such effective models can be applied even for densities where 
the underlying theory, QCD, becomes to leading order a free Fermi gas. 

Back to the stability analysis: Note that the signs of the static eigenvalues g 
therefore the stability of the uniform regime follow from the signs of the derivatives 
acting on Q and from the signs of the coefficients Xi^N,o (H)- From all possible Xi^N,a 
with z = 1, 2, 3, > 1 and a = 0, 1 only the term Aj=i_Ar=i^a=o is negative, all the 
other coefficients are guaranteed to be larger than zero or at most equal zero. (To the 
latter category belong the six 0(4) zero modes with N = 1 and a = 1. Furthermore 
all the terms Aj=3^Ar^i are zero indicating that the sixth-order term itself has infinitely 
many zero modes.) Now taking the high-density (Q) and the low-density behavior 
(|3^) into account, we can conclude that at high densities the uniform B = 1 skyrmion 
solution on the hypersphere is bound to be stable whereas at sufficiently low densities, 
where the non-linear a model term becomes the dominant one, the uniform solution 
is unstable since A^^^^ q becomes negative eventually. Furthermore we have to take 
into account that in the stable regime - even infinitesimally close to the instability - 
all A^ have to increase monotonically with N, that the A^ i modes have to be stable 
by fiat and that the eigenvalues of the modes have to be smooth functions of L. Then 
the existence of a critical value of L is guaranteed where A^ g becomes negative such 
that the uniform high density phase becomes unstable. There is a bifurcation from 
the uniform solution which is chirally restored to the usual, localized and chirally 
broken B = 1 hedgehog solution. So even for the most general geometric lagrangian 
(|2^) built form first order field derivatives the existence of the chiral phase transition 
(as discussed at length in the preceeding chapters for the normal Skyrme model) 
is guaranteed under the assumption of a few reasonable constraints (the vacuum 
structure, Eqs. ( ^3]) and (0)) on the form of the lagrangian. 



4.2. Vector Meson Stabilization 

We have discussed so far lagrangians which ensured the stabilization of skyrmions 
by higher-order (contact) terms in the field derivatives. Naturally the question arises 
whether the existence of the chiral phase transition is linked to the stabilization by 
contact terms or whether vector mesons can be added which act over a finite range, 
as e.g. uj (and p) mesons. Let us take as an example a Skyrme- type variant with 
uj meson stabilization. (See ref.li^ for a periodic array study of this model.) The 



corresponding lagrangian has the structure Wi. 

C2,u. = ^Tr(9^[/t9^f/) - -u^.u^" + -fu^u^ + g^u^B^- (35) 

where the first part is the usual non-linear a model term, the second one the uj kinetic 
term expressed through the u) field-strength tensor uj^y = d^ujy — d^u^, the third part 
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is the uj mass term and the last one is the couphng of the uj field to the topological 
baryon current As shown by Adami El the lagrangian (^) has stable B = 1 

solitons. 

By putting an cj-stabilized B = 1 hedgehog on the hypersphere one can easily find 
that one solution of the equations of motion is always the uniform solution /(//)= 
for the hedgehog profile and 

for the cuq- component. The Wo-component is in this case spatially constant and de- 
pends only on the hypersphere radius L. (As in fiat space the spatial components 
uji are identically zero because of the static hedgehog form of the soliton profile.) 
When the uniform solutions /(/i) = /i and (^) are reinserted into the corresponding 
energy density of the lagrangian (pS]), one recovers the same structure as for a model 
stabilized by a sixth-order term with = —{g'^/2mlj)B^B^. This should not come 
as a surprise, since in the limit m^,gi^ —>■ oo, under a constant ratio Quj/^uji the la- 
grangian (P^D converges to the lagrangian (123) with cg = Both lagrangians 



have the same high-density phase in the static sector. The physical reason for the 
identical behavior of the static solutions of both lagrangians is that at high densi- 
ties it is energetically favorable that the tuo-component becomes spatially constant, 
otherwise the "costs" in energy from a finite spatial gradient term in this component 
becomes higher and higher with increasing density. This together with the equation 
of motion for the cuo-component necessarily lead to Eq.(^) which in turn guarantees 
the identical high-density behavior of both lagrangians. However, the cj-stabilized 
model bifurcates at a smaller and the low-density behavior of both lagrangians is 
different. Under increasing parameters g^^ and with the ratio gui/ni^^ kept constant 
the critical density and the low density phase of the uj stabilized model will approach 



the corresponding quantities of the model (24). So the inclusion of finite range vector 



mesons leads to an increase in the transition density in comparison to the correspond- 
ing contact-term model, see also ref.O. In summary, also the cu-stabilized Skyrme 
model has a critical hypersphere radius where the uniform (chirally restored) phase 
on the hypersphere becomes unstable and bifurcates into a localized (chirally broken) 
phase which has the usual B = 1 hedgehog as fiat space hmit. 

A similar behavior follows when the £2,lj niodel is extended by the introduction 
of p mesons. This can be done in different ways. For our purposes the only essen- 
tial precondition is that the total lagrangian is still explicitly chiral symmetric. For 
si mpli city let us consider in the following the hidden p meson coupling a la Bando et 
al. in a minimal way 1^, i.e. 

C2,u,p = ^Tr(9^ma^f/)-a^Tr(eta^e-e9^e^-2z^7P^)2-ip^,p^^ 

1 TTL^ 

- -uj^yUj^" + -^u^uj^ + gujUJfj.B'' + {ujpiT coupling terms) (37) 
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with ^ = v^, and p^^, = d^j^py — d^p^ — ig[p^,pj\, the non-abehan field-strength 
tensor. The second term in ( |57D is responsible for the generation of the p mass and 
the pTTTT coupling with a coupling constant g. (The standard choice for the parameter 
a is a = 2, see ref. El.) Note that under the static hedgehog ansatz only the spatial 
Pi components are excited. Again it is easy to show that the equations of motion of 
this model are satisfied by the uniform profiles 

/(/i) = p (38) 
uo = -^Bo (39) 

A = with e = exp(^r-f/x/2). (40) 



The last equation guarantees that the second term in ( p7D vanishes and that the p 
kinetic term behaves as fourth-order Skyrme term with a coefficient el = 1/ [Sg"^). The 
lagrangian ( p7| ) (without extra upir couplings) has therefore the same high density 
behavior as the simple contact lagrangian 



A,4,6 = ^TT{d,Wd^U) + -L_Tr[Wd,U,WdM? - I (y^Y B^B^. (41) 



This behavior presented here for the simplest djpvr-model is generic - just the co- 
efficient of the sixth-order term in (^T]) may change when tupvT-coupling terms are 
introduced. For all models with vector-meson stabilization which have been studied 
far the following is true: It can be numerically shown that 



so 



(i) the only existing B = 1 solution at high densities is the uniform (delocalized) 
hedgehog one, the static solutions are identical to those of contact-term la- 
grangians with suitable coefficients and the static vector meson fields are equal 
to their driving pionic currents {see e.g. (^) or (plj)), 

(ii) there exists a critical model-dependent hypersphere-radius where the uniform 
solution becomes unstable and bifurcates into the usual localized configuration 
and 

(iii) at low densities the solutions on the hypersphere S^(L) converge asymptotically 
to the corresponding solutions in the fiat space. 

The essential conditions on these models are two-fold: they should be explicitly chiral 
symmetric and they should allow for locally stable B = 1 hedgehog solitons. In 
case these two conditions are met, there will be necessarily a chiral phase transition 
at high densities with the same features as already discussed. Still, under realistic 
values for the mesonic input parameters {e.g. /,r, £4, cq, see ref. IH, or (7^ and g, see 
refs. S H' i)) the transition densities pc come out systematically too low, pc ~ Pnm 



(0.16 fm~^). Apparently these models lack terms which generate additional attraction. 
Note that the hedgehog mass under such realistic parameters is systematically too 
high (~ 1.6 GeV) whereas there is not sufficient central attraction in the B = 2 



16 



system.! If one would tune the stabilizing parameters - while keeping = 93MeV 
fixed - such that the classical hedgehog mass would be 0.87 GeV (consistent with a 
nucleon mass of 0.94 GeV), then the transition density pc would be approximately 
three times the nuclear matter density. 

5. Second- Versus First-Order 

Whereas in the periodic array calculations second as well as first order phase 
transitions were found depending on the choice of the lattice and the twisted boundary 
conditions, we have so far encountered only second order transitions between the 
delocalized and the localized phase on the 3-sphere. This does not need to be the 
case in general: All the Skyrme-type lagrangians presented are explicitly not scale 
invariant and do not "know" (yet) about the trace anomaly in QCD. Let us therefore 
try to incorporate the same scaling behavior as in QCD into these lagrangians (see 



refs. 1^ ^ B3), e.g. consider the lagrangian 

A,4,x = Tr(a^f/t9'^t/) + |Tr[[/t5^[/,[/t9.t/]2 + ^(9^x)' 




-5b 1+ - log (42) 



as simplest extension of the usual Skyrme lagrangian. The scalar field x with the 
vacuum expectation value xo = (0|x|0) is introduced with the purpose of making the 
first term in (^) scale invariant (the Skyrme term and the x-kinetic term are already 
scale invariant), whereas the last term in (^), the x-potential term, is adjusted to 
fit the trace anomaly of QCD. Bb is the "bag constant" which can be expressed 
in terms of the gluon condensate as Bb = (9/32)(0|(as/7r)G'^|0). For values of the 
parameters Xo and Bb see ref. 0. (Note that the fluctuations of the x-field correspond 
to glueball-excitations.) In ref. EH it was numerically shown that the lagrangian (^) 
again possesses a phase transition to a uniform delocalized solution at high densities 
with /(/i) = /i and a constant x profile which - depending on the choice of parameters 
- can be even x = 0- Because of the insufficient accuracy in the numerics, the authors 
of ref. £3 missed the fact that the phase transition is of first order. The latter is more 
or less obvious from the form of the x-potential term which is adjusted to be minimal 
for X = Xo- Thus it is impossible to find a smooth transition between the constant 
X < Xo (high-density) profile and the localized (low density) x profile (which at the 
south pole on the 3-sphere is exactly xo) without violating the x equation of motion: 



sm jjL Xo V A* / 



^11 V I V \ 

(43) 




For a non- vanishing nonlinear sigma term (as in the present case since /(p) 7^ 0) 
the X field cannot both be constant and equal to xo- So the x field can only do the 
transition from the low- to the high-density phase and vice versa by a jump. Thus by 
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incorporating a scalar field a la Schechter H in the Skyrme model (or its extensions) 
the second order phase transition can be easily changed to a first order one. See also 
ref. E3 for a different mechanism to achieve the same thing. 

It is rather easy to extend the lagrangian ( P?]) to the most general with the QCD 
trace anomaly consistent geometrical form involving only first derivatives to even 
order: 




+ lid.xr -Bb\1+{^] log ( ) I . (44) 
Note that the restoration of scale invariance, x = 0, is only consistent with the 



lagrangian (0) if Q has the asymptotics (0), i.e. 

I G (-i-f^2, -(-)'A, ^ Ci-C^y + C2(-A) if X - 0. 

\ X X X J 

(45) 

We know from what was said before that one can expect in this case a first-order phase 
transition to a chirally restored phase. However, the x-profile in the high-density 
phase does not need to become zero, it can equally well be just a constant x < Xo- 
Scale invariance implies chiral restoration, but not vice versa. Nevertheless, just the 
possibility to have a scale invariant limit at high densities leads to strong constraints 
on the form of the effective model. If a lagrangian is considered which includes any 
term higher than fourth-order in the derivatives of the pion fields, terms of all orders 
have to be included as well both to satisfy the high-density asymptotic behavior (^) 
and to ensure the stability of the chiral symmetric high-density phasclH The necessity 
for infinitely many terms is consistent with the large N^. philosophy.!!' H Because of 
asymptotic freedom, the effective model which follows from QCD to leading order in 
the 1/Nc expansion must involve infinitely many mesons. The hope is - as shown 
for the u and p meson case - that the high density phase of this infinite tower of 
meson resonances can be approximated by infinitely many suitable contact terms. 
This might lead to "dreams" of Regge trajectories. 

6. Discussion 

Even effective models of the Skyrme class (which do not possess any quark degrees 
of freedom) can have a phase transitions to a high-density phase with restored chiral 
symmetry. There are essentially just two conditions on such a model: It must have the 
relevant symmetries (especially the lagrangian must be explicitly chirally invariant) 
and it must treat the baryon structure and interaction on the same footing. Normally, 
the phase transitions are of second order, but one can easily change them to first- 
order ones by incorporating the QCD trace anomaly for instance. This allows also - 
depending on the input parameters - for a scale invariant high density limit. Scale 
invariance implies chiral symmetry restoration, but not vice versa. 
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Under the same parameter input (from the meson sector) the periodic array cal- 
culations as well as the 3-sphere calculations give approximately the same values for 
the transition densities. The predictions, however, are for both systematically too 
low, pc ~ Pnm- This problem may be connected with two other deficiencies which 
plague the Skyrme-like models: the missing central attraction in the B — 2 channel 
and the predicted high value for the hedgehog mass. Quantum (higher loop or higher 
1/Nc) corrections may have a chance to cure all three problems at the same time by 
providing extra non-local terms which generate attraction. 
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